Cooling the Collective Motion of Trapped Ions to Initialize a Quantum RegisterQ 
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We report preparation in the ground state of collective modes of motion of two trapped 9 Be + 
ions. This is a crucial step towards realizing quantum logic gates which can entangle the ions' 
internal electronic states. We find that heating of the modes of relative ion motion is substantially 
suppressed relative to that of the center-of-mass modes, suggesting the importance of these modes 
in future experiments. 
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Q\ • In physics, quantum computation Jy provides a general framework for fundamental investigations into subjects 
ON , such as entanglement, quantum measurement, and quantum information theory. Since quantum computation relies 
on entanglement between qubits, any implementation of a quantum computer must offer isolation from the effects 
(3J[)! of decoherence, but also allow controllable and coherent interaction between the qubits. Cirac and Zoller || have 
^ . proposed an attractive scheme for realizing a quantum computer, which is scalable to an arbitrary number of qubits. 
Their scheme is based on a collection of trapped atomic ions, where each qubit (one per ion) is comprised of a pair of 
the ions' internal states, while quantum information is transferred between different ions using a particular quantized 
mode of the ions' collective motion. This "quantum data bus" must first be initialized in a pure quantum state [p|: 
for example, its ground state The basics of this scheme have been demonstrated experimentally in a fundamental 
t^J- 1 logic gate (a Controlled-NOT) operating between a motional mode of a single trapped ion and two of the ion's internal 
\ states ||. In that work, the motional state was initialized in the ground state by laser cooling ||. The next step 
C*") i towards implementing the Cirac-Zoller scheme is to cool at least one mode of collective motion of multiple ions to the 
ground state. In this Letter, we describe the first experiments to realize this goal. We also report significant difference 
, between the decoherence rates of the center-of-mass and non-center-of-mass modes of motion. 

We confine 9 Be + ions in a coaxial- resonator-based rf (Paul) trap, similar to that described in Ref. ||. The electrodes 
00 , in this trap are made from 125 /um-thick sheets of Be metal, as shown in Fig. [j]. We apply a potential <fi{t) = 
Vq cos(r2x^) + Uq to the (elliptical) ring electrode relative to the endcap electrodes. If several ions are trapped and 
cooled, they will naturally align themselves along the major axis of the ring electrode. The electrode's elliptical shape, 
in combination with Uq > 0, allows a linear crystal to be maintained while suppressing rf-micromotion of the ions along 
this direction JjJ. With Vq « 520 V, Qt/2tt ss 238 MHz, and Uq —0 V, the pseudopotential oscillation frequencies are 
(w x ,Wy,w z )/2Tr w (4.6,12.7,17.0) MHz. With Uq = 18.2 V, the frequencies become (8.6,17.6,9.3) MHz. Fig. | shows 
two ions confined in the trap and imaged with an //3 lens system onto a position-sensitive photomultiplier tube. 
The ions are cooled and probed with laser beams whose geometry is indicated in Fig. |](a). The relevant level 
O h> structure of 9 Be + is shown in Fig. ||(b). The quantization axis is defined by an applied static magnetic field; \B\ « 
0.2 mT. The levels of interest for quantum logic operations are the 2s 2 Si/ 2 \F — 2,mp — 2) and 2s 2 Si/ 2 \F — 1, ttif = 
1) states, abbreviated by | |) and | f), respectively. Laser beams Dl, D2, and D3 are cr + -polarized and focussed to 
nearly saturate the ions (I sa t ~ 85 mW cm -2 ). Beams Dl and D2 provide Doppler precooling in all three dimensions, 
and beam D3 prevents optical pumping to the \F = 2,mp = 1) state. The | |) — > 2p 2 P 3 / 2 |F = 3,itif = 3) transition 
(radiative linewidth 7/27r 19.4 MHz), driven by D2, is a cycling transition, which allows us to detect the ion's 
electronic state (| J.) or | f)) with nearly unit detection efficiency. 

Beams Rl (o^/er^-polarized) and R2 (7r-polarized) are used to drive stimulated Raman transitions between | J.) 
and | |), through the virtual 2p 2 P 1 / 2 state p). These beams are derived from a single laser, whose output is split by 
an acousto-optic modulator ||. The beams are detuned by A/2tt w 40 GHz to the red of the 2s 2 S , 1 / 2 —* 2p 2 Pi/2 
transition, and their frequency difference is tuned around the 2s 2 Si/ 2 hyperfine splitting of luq/2it ss 1.25 GHz. 
(Here, ojq includes stable shifts of a few megahertz from the Zeeman and ac Stark effects.) R2 is directed along 
(-1/V2)x + (l/2)(-y + z). If R1_LR2 as in Fig. ||, then the Raman beam wavevector difference 8k || x, and 
the transitions are sensitive to ion motion only in this direction. If, however, Rl is counterpropagating to R2, the 
transitions become sensitive to motion in all three dimensions. 

When two cold ions arc held in the trap and undergo small oscillations about their equilibrium positions, we may 
solve the equations of motion using normal mode coordinates. For two ions lying along the x-axis there are two modes 
involving motion along this axis: the center-of-mass (COM) mode (in which the ions move together with frequency 
wcom = u x ) and the stretch mode (wherein the ions move out of phase, with frequency w str = V^^com)- The other 

motional frequencies are oj v (y center-of-mass), ui z (z center-of-mass), \/uj 2 ~lu 2 (xy rocking), and uj 2 — u> 2 (xz 
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rocking). 

The lower traces in Fig. |^, taken with 5k \\ x, shows an x-axis normal mode spectrum; results for the y- and z-modes 
are very similar. We take the data with the following steps: first we turn on beams Dl, D2, and D3 for approximately 
10 /is to Doppler cool the ions to the Lamb-Dicke regime, where the ions' confinement is much smaller than the laser 
wavelength. Next, we turn off beam D2, and leave beams Dl and D3 on for 5 /is to optically pump both ions to 
the | [) state. We then turn on only the Raman beams Rl and R2 for a time t pr , with relative detuning luq + 5 pr 
(the "Raman probe" pulse). Finally, we drive the cycling transition with D2 and measure the ions' fluorescence. We 
repeat the experiment at a rate of a few kilohertz while slowly sweeping 5 pr . If the Raman beam difference frequency 
is resonant with a transition, then an ion is driven from | 4) — ► | T) and the D2-driven fluorescence rate drops. 

For a single ion, the carrier transition (5 pr — 0) causes the population to undergo sinusoidal Rabi oscillations 
between | 4) and | 4) ||. The effective Rabi frequency is Q — gi^/A ~ 2ttx 250 kHz, where g\, g2 are the single- 
photon resonant Rabi frequencies of beams Rl and R2. (We assume A 3> 7, u) m 3> where u> m is the frequency of 
the motional mode of interest.) If 5 pr = —uj x (the first lower x sideband), then the transition couples the states | 4, n x ) 
and | ~\,n x — 1), where n x is the vibrational level of the quantized motion along x. In the Lamb-Dicke regime, the 
corresponding Rabi frequency is given by r2n x ,n x -i = rjxy/n^tl ||. Here, r\ x — xo\Sk • x\ is the Lamb-Dicke parameter 
(= 0.23 when uj x /2ir = 8.6 MHz) and Xq = yhj (2muj x ) is the spread of the n x = wave function, with m being the 
ion's mass). (Note that if the ion is in the n x = state of motion, this lower sideband vanishes.) The first upper x 
sideband transition (S pr — +w x ) couples | i,n x ) and | 1,n x + 1) with Rabi frequency ftn^.n^+i — , Ux\Jn x + lf2. 

In the case of two ions driven on the carrier transition, each ion independently undergoes Rabi oscillations between 
| 4) and | 4) with Rabi frequency f2. Since the laser beam waists (~ 20 /im) are much larger than the ion-ion 
separation (w 2/im), the ions are equally illuminated. Nonetheless, if the micromotion of the two ions is different, 
then the reduction of the carrier (and sideband) transition strengths due to the micromotion will giv e a different Rabi 
frequency for each ion J? 10 1. This could be used as a means of selectively addressing the ions 11 1; however, in the 
present work the two ions' Rabi frequencies were equal. 

Since the sideband transitions affect the motional state, which is a shared property of both ions, such transitions 
produce entanglement between the ions' spins and their collective motion fl2| . The system can no longer be treated 
as two, independent, two-level systems and the measured fluorescence following the Raman probe is a complicated 
function of the probe pulse duration t pr . For example, given an initial state | 4, 1, n) (where n is the vibrational level 
of the COM or stretch motion along the x-axis) driven on the corresponding lower sideband for a time t pr , the wave 
function evolves as 

r n ^ I n (\ l>t)±e # | ti))\n-l) 
\Mt P r)) = [l - ^— y [1 - cos(GV)] ] 14, In) - ^ (e ^ )/2 \I^Y sin(GVp -—L 

^ ^n-l ^ 1 ~ cos ( G V)] I 4, 4, n - 2), (1) 



where G = W2{2n — 1)0 r\ x . m and 9, cf> are the sum and difference of the Raman beam phases at the ions. On the 
COM sideband (top sign in Eq. (|lf)), r\ x . m = i?i.com = Vx/V2 (down by a factor of \[2 from the single-ion case due 

to the extra mass of the two-ion string), whereas on the stretch sideband (lower sign), r\ x , m = ^str = Vx/v 2v^3- 
The expressions for transitions on the upper motional sidebands are similar. If, before the Raman probe pulse, the 
ions have probability p n of being in the motional state |n), the subsequently- measured average fluorescence from the 
cycling transition is 

S(t pr ) = ]>>n (2|(4,4,n|^„(v))l 2 + 1(4, T,n- i|Vn(V))l 2 + l(T,4,n- i|^„(V)>l 2 ) ■ (2) 

n 

This signal is proportional to the expectation value of the number of atoms in the state | 4)- For the data shown in 
Fig. ||, t pr was chosen to maximize the sideband features. 

The upper traces in Fig. || show the effects of adding several cycles of Raman cooling || on one particular x-mode 
after the Doppler cooling but before the probe pulse. The reduction in the mean vibrational number (n) is indicated 
by the reduction in size of the lower sideband, which vanishes in the limit (ri) — > 0. The data are consistent with 
a thermal state of (ticom) = O-lll^Jg or (n stT ) — O.OIIq'q*. This implies that the COM and stretch modes are in 
their ground states 904^2% and 99^% of the time, respectively. We have also simultaneously cooled the COM and 
stretch modes along x, to comparable values of (n) (and have separately cooled the other four motional modes — y 
and z COM, xy rocking, and xz rocking — to near their ground states). 

Each cycle of Raman cooling consists of: (1) a pulse of the Raman beams with their difference frequency tuned to 
one of the lower sidebands (COM or stretch mode) and (2) optical repumping to the | 4) state driven by beams Dl 
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and D3. The Raman transition reduces the vibrational energy by huj m , whereas the repumping, on average, heats 
each mode by approximately the recoil energy (<C fiw m ). Therefore, the ion is cooled through this process. Five pulses 
of Raman cooling were used for the data shown in Fig. || The exact durations of the Raman pulses were chosen to 
optimize the cooling rate each pulse was approximately 5 [is long. 

For an ion-trap implementation of a quantum computer, the motional modes are most susceptible to decoherence. 
The ions' motional states lose coherence if they couple to (stochastic) electric fields caused by fluctuating potentials 
on the electrodes. This leads to heating, which has previously been observed in single ions |[5 
heating drove the ion out of the motional (COM) ground state in approximately 1 ms. We have performed similar 
heating measurements on the COM and non-COM modes of motion of two ions. The results are summarized in 
Table |. The heating rate was determined by inserting a delay between laser cooling and the Raman probe. The 
main results from these data are that the COM modes are heated out of the ground state much more quickly than 
the non-COM modes. This can be explained as follows. 

The COM modes, in which both ions move in phase, can be excited by a uniform electric field. However, no 
non-COM mode can be excited by a uniform electric field — since these modes involve differential motion of the 
ions, they can only be driven by field gradients. If the fluctuating field at the ion (along the direction of motion of 
the mode of interest) is E(t), an estimate of the corresponding field gradient is E(t)/d, where d is a characteristic 
internal dimension of the trap. For stochastic fields, the COM heating rate scales as (E 2 (t)); the non-COM mode 
heating rates scale as ([^jp Ax] 2 ) (where Ax is the ion-ion separation), down by a factor of I0 4 for the present trap. 
Similarly, other non-COM modes for more than two ions can only be excited by higher-order field gradients, leading 
to further reductions in their heating. 

This suggests using non-COM modes for the quantum data bus in the Cirac-Zoller scheme. Excitation of the 
"spectator" COM modes along the direction to which the Raman transitions are sensitive will still alter the Rabi 
frequencies, but these effects will be higher order in the Lamb-Dicke parameter In the two-ion example, in 

the Lamb-Dicke regime, the Rabi frequency for a first sideband transition \m) — ► \n\) on (cold) mode 1, given that 
(hotter) mode 2 is in the state \n 2 ), is Q 

(n 2 ) = n m ^Te- {r ' 2 ^ )/2 (l - n 2 r,l), (3) 

where n 1> denotes the larger of or n\, and 771 and 772 are the Lamb-Dicke parameters for modes 1 and 2, respectively. 
Fluctuations in the Rabi frequency of mode 1 due to fluctuations in n 2 therefore occur in order r] 2 . However, for the 
conditions of the present experiment, even if quantum logic operations were performed using the x-stretch mode, the 
i-COM mode heating would still limit the number of operations to around ten by the above mechanism. Clearly, this 
heating must be eliminated in future experiments. 

The two- ion cooling results presented here are comparable to our previous single-ion results [[|, indicating that rf 
heating should not be a concern for small numbers of ions Jic| l . Comparable cooling for N > 2 ions should not present 
any fundamental difficulties, as long as spurious overlaps of motional modes are avoided. 

The preparation of a pure state of motion (the ground state) of multiple trapped ions represents the first step towards 
realizing quantum logic operations on them. Such operations should lead to the creation of arbitrary entangled states 
of massive particles, including EPR- or GHZ-like spin states Unlike other systems in which EPR states have 

been generated, it should be possible to reliably create these states on demand [[ll] rather than by a selection process, 
and to detect them with nearly perfect efficiency p6[ . 
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FIG. 1. Two ions trapped in an elliptical rf (Paul) trap. 
The ring has an aspect ratio of 3:2 and the major axis is 
525 /im long. The slot which forms the endcaps is 250 fim 
across. A potential cj>(t) is applied to the ring (see text). The 
Be sheets are ~ 125/im thick. With an :r-axis pseudopotential 
oscillation frequency lu x /2it — 4.6 MHz, the ion-ion spacing 
is approximately 3 /im. 




FIG. 2. (a) Laser beam geometry. The trap ring elec- 
trode is shown rotated 45° out of the page. The endcap elec- 
trodes are omitted for clarity; see Fig. hi A magnetic field 
B of magnitude 0.2 mT defines the quantization axis along 
-{l/y/2)x + (1/2) (tf - z), and laser beam polarizations are 
indicated, (b) Relevant 9 Be + energy levels (not to scale), in- 
dicated by F, mf quantum numbers in the ground state. 2 P 
fine structure splitting is « 197 GHz, 2 Si/2 hyperfine splitting 
is lo /2tt » 1.25 GHz, 2P 1/2 hyperfine splitting is « 237 MHz, 
and the 2 P3/2 hyperfine structure (<C 7/27T « 19.4 MHz) is 
not resolved. All optical transitions are near A ~ 313 nm, and 
A/2?r « 40 GHz. 



FIG. 3. Spectrum of sidebands due to two-ion a;-axis nor- 
mal mode motion: (from left to right) lower stretch, lower 
COM, upper COM, and upper stretch. The ordinate is the 
detuning of the Raman probe beam difference frequency from 
the carrier transition. The abscissa shows the ion fluorescence 
(proportional to the expectation value of the number of atoms 
in the state | J,)), plus a constant background (whose approx- 
imate level for the lower curves is indicated by the dashed 
line). The solid lines, meant as guides to the eye, are fits 
to Gaussians. The lower traces show the effects of Doppler 
cooling. The upper traces, offset vertically for clarity, show 
the effects of several pulses of Raman cooling on the mode 
which is displayed. Vanishing lower motional sidebands indi- 
cate cooling to the ground state of motion. The peak widths 
are consistent with the Raman probe pulse lengths (~ 3 fis). 

TABLE I. Heating rates of the six normal modes of two 
trapped ions. The Raman beams were counterpropagating for 
the y- and z- axis data, making the Raman probe sensitive to 
motion in all three dimensions. Note that the COM modes 
are heated at a much higher rate than the non-COM modes 
(see text). (The precision with which the heating rates are 
given for the last five modes is limited by measurement noise.) 
mode w m /27r (MHz) 5{n)/5t (ms _1 ) 



xcom 8.6 19t$ 

2/com 17.6 > 10 

zoom 9.3 > 20 

X s tr 14.9 < 0.18 

Crocking 15.4 < 1 
% ^rocking 

3.6 < 0.5 
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